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Application of the fractional calculus to quantum processes is presented. In particular, the quan- 
tum dynamics is considered in the framework of the fractional time Schrodinger equation (SE), 
which differs from the standard SE by the fractional time derivative: ^ -g^. It is shown that for 
a — 1/2 the fractional SE is isospectral to a comb model. An analytical expression for the Green 
functions of the systems are obtained. The semiclassical limit is discussed. 
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Application of the fractional calculus to quantum 
processes is a new and fast developing part of quan- 
tum physics which studies nonlocal quantum phenomena 
[U, H i, i, i, S 0, i, i, ■ It aims to explore non-local 
effects found for either long-range interactions or time- 
dependent processes with many scales [1, [ll|, [H, [H, [13, 
[l5, 16]. Fractional calculus with a variety of applications 
H, [12|, III, [13, III [li, [ll M^Mb including amlications 
to quantum processes 0, UHUi H [3 0, HH, [H, HJ] , 
is a well developed and well established field that was 
extensively reviewed. 

The concept of differentiation of non-integer orders 
rises from works of Leibniz, Liouville, Riemann, Grun- 
wald and Letnikov, see e.g., [H, U^. Its application is 
related to random processes with power law distributions. 
This corresponds to the absence of characteristic average 
values for processes exhibiting many scales [U, [2^. A 
continuously increasing list of applications in many sci- 
ences has developed: it includes material science y 
physical kinetics f\A] , anomalous transport theory with a 
variety of applications in solid state physics [111, [H, [25j 
and in nonlinear dynamics [3, ■ 

In quantum physics, the fractional concept can be in- 
troduced by means of the Feynman propagator for non- 
relativistic quantum mechanics as for Brownian path in- 
tegrals [131. Equivalence between the Wiener and the 
Feynman path integrals, established by Kac [1^, indi- 
cates some relation between the classical diffusion equa- 
tion and the Schrodinger equation. Therefore, an appear- 
ance of the space fractional derivatives in the Schrodinger 
equation is natural, since both the standard Schrodinger 
equation and the space fractional one obey the Markov 
process. As shown in the seminal papers [2, [3]' it relates 
with the path integrals approach. As a result of this, the 
path integral approach for Levy stable processes, leading 
to the fractional diffusion equation, can be extended to 
a quantum Feynman-Levy measure which leads to the 
space fractional Schrodinger equation [3, [3] • 

A fractional time derivative can be introduced in 
the quantum mechanics by analogy with the fractional 
Fokker-Planck equation (FFPE), as well, by means of the 
Wick rotation of time t —it/h |8|. However its physical 
interpretation is still vague: for example, a phase of the 
wave function as well as the semiclassical approximation 
should be understood. The fractional time Schrodinger 



equation was first considered in [8| . Its generalization to 
space-time fractional quantum dynamics [13j[lll was per- 
formed and a relation to the fractional uncertainty (2^ 
was studied as well. Exact solutions were also obtained 
for the time fractional nonlinear Schrodinger equation 
[23j . It is worth noting that, contrary to the space frac- 
tional derivative, the fractional time Schrodinger equa- 
tion describes non-Markovian evolution with a memory 
effect. 

The fractional time quantum dynamics with the 
Hamiltonian H{x) is described by the fractional 
Schrodinger equation (FSE) 



iJV(x, t) 



(1) 



where a < 1. For concordance of the dimension in Eq. 
P]) all variables and parameters are considered dimen- 
sionless, and h is the dimensionless Planck constant, see 
also [1, 20]. For a = 1, Eq. ^ is the "conventional" 
(standard) Schrodinger equation. For a < 1 the frac- 
tional derivative is a formal notation of an integral with 
a power law memory kernel of the form [29j 



dt 



r(i 



dr 



dr, (2) 



which is known as the Caputo fractional derivative [r7| . 
and r(z) is a gamma function. This definition makes 
it possible to carry out the Laplace transform of the 
fractional derivative. Introducing the Laplace image 
4'{s) = C4>{t), one obtains 



C 



V(0), 



(3) 



Another interesting property of the FSE is time evo- 
lution in the form of the Mittag-Leffler function. For 
the time-independent Hamiltonian, the eigenvalue equa- 
tion (with corresponding boundary conditions) is iJ(/)A = 
A</>A. Therefore one obtains the Green function in the 
term of the Mittag-Leffler function 1(2) = Ea{z) = 
E^o^Vr(ja-t-l) [i,[il: 



G(x,t;x')=^<^I(x')0A(x)i?„ 
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which is a fractional generaHzation of Green's function. It 
is worth noting that this solution for the Green function 
in the form of the Mittag-LefHer function does not satisfy 
Stone's theorem on one-parameter unitary groups [scf]. 

In the general case, when the eigenvalue problem can- 
not be solved rigorously, the analysis of the fractional 
Green function meets serious deficiencies. For a exam- 
ple, the semiclassical analysis of the FSE leads to a much 
more complicated form for the wave function than the 
physically transparent expression for the local wave func- 
tion il){x,t) ~ giS{t,x)/h^ Moreover, the classical action 
S{t, x) is not defined anymore for the FSE. Another im- 
portant question could be addressed to both FSE ([T]) and 
Eq. (jlj: for the fractional Fokker-Planck equation, the 
fractional time derivative is an asymptotic description of 
a continuous time random walk and it describes subdif- 
fusion. Therefore, what is the physical meaning of the 
fractional time derivative in quantum processes when a 
concept of a multi-scale continuous time random walk is 
absent? 

To shed light on this situation, we consider a case with 
a = 1/2, when fractional quantum dynamics can be mod- 
elled by means of the conventional quantum mechanics 
in the framework of a comb model. The comb model 
is an analogue of a Id medium where fractional diffusion 
has been observed [Sll, [12] ■ It is a particular example of a 
non-Markovian phenomenon, explained in the framework 
of a so-called continuous time random walk 
This model is also known as a toy model for a porous 
medium used for exploration of low dimensional percola- 
tion clusters 

A special quantum behavior of a particle on the comb 
is the quantum motion in the d -I- 1 configuration space 
(x, y) , such that the dynamics in the d dimensional con- 
figuration space X is possible only at j/ = and motions 
in the x and y directions commute. Therefore the quan- 
tum dynamics is described by the following Schrodinger 
equation on a comb 



(5) 



where the Hamiltonian is the same as in Eq. ([T]), and 

H = H{x) = + V(x) governs the dynamics with 
a potential V^(x) in the x space, while the y coordi- 
nate corresponds to the Id free motion. All the pa- 
rameters and variables are dimensionless [s^. We will 
study an initial value problem with the initial condition 
^(t = 0) = 5'o(x, y). It is worth noting that the semi- 
classical asymptotic expansion in h of the wave function 
is impossible, since the S potential in the y direction, 
"cannot have a sensible semiclassical limit" [35, 36] and 
degrees of freedom cannot easily be separated. Using the 
eigenvalue problem 



(6) 



we present the wave function in Eq. ^ as the expansion 
'i'{x,y,t) = X^A '/'A V'a(x), where J2x also supposes 



integration on A for the continuous spectrum. For the 
fixed A we arrive at the dynamics of a particle in the S 
potential 



ih- 



dt 
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The Green function for this Schrodinger equation has 
been obtained in [H, \3T\ 

Gx{y,t;y')^Go{y,t;y') + 



duGo{\y\ + \y'\ + u,t:0) exp(-MA//i)) , (8) 
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where Gofy, t; y') = . ^ exp \ i j ^ ) is the free par- 
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tide propagator. Using this result, we obtain for the wave 
function of Eq. (O 

*(x, y, t) = / dy'Goiy, 2/')*o(x, y') + 
- J dy'duGo{\y\ + \y'\+ u, t- 0)e-"^/''f vI/o(x, y') . (9) 

Now our aim is to compare Green's function of this so- 
lution with the one of Eq. @ . Contrary to fractional dif- 
fusion, where the FFPE with a = 1/2 is identical to the 
comb model, the FSE ^ is not identical to the quantum 
comb model of Eq. ([5]). Nevertheless, the models have 
some features in common, namely there are the some 
singularities of the Green functions which determine the 
spectrum. To show this, let us present both Eq. ^ and 
Eq. ^ in the form of the inverse Laplace transform. 
For simplicity we consider the one-dimensional x space. 
First, we consider Eq. ([9]). Presenting the wave function 
as the Laplace inversion and carrying out integration on 
u in the second part, we have for the Green function 

G{x, y, t; x', y') = Go(y, t, y')5{x - x') - C-Z Ea 

oxp (^st+yj-2is/h\\y\ + \y'\\jds 



-**(x')«'a(x), (10) 



where the eigenvalue problem of Eq. ([6]) is used. Using 
integral presentation of the Mittag-LefHer function [isl ]. 
we rewrite the Green function of the FSE in Eq. ([4]) as 
follows 



cr+ioo 



(11) 



(T-ioo ^[y/s — X/ V2ih) 



Therefore, disregarding the free propagator in Eq. ([Tt 
one obtains the same spectral properties of both systems, 
since the spectral decompositions of the evolution oper- 
ators are the same. 

In what follows we consider the dynamics of the FSE in 
the framework of the Green function (jlip . There are two 
contributions to this integral from a pole at sg — X^ /2ih 
and from a branch point at s = with a branch cut from 
— oo. Carrying out this integration (see @), we present 
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the wave function in the following operator form (see Eq. 
®) 



tively) is nonlinear and the expectation values can di- 
verge for the same initial condition i/jqIx): 
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It consists of two parts, the oscillatory one and the decay 
in time. This expression is convenient for studying the 
large time asymptotic t ^ 1, when the decay term can 
be neglected. 

The decay term is of the order of ^ V^, and can be 
neglected in the semiclassical limit ft, as well. The 
evolution of the semiclassical wave function ^sc\{x^t) = 
tfj{x,t)/2 is due to the oscillatory term with the Hamil- 
tonian Ti,sc\ = This semiclassical dynamics is es- 

sentially nonlinear. 

As an example of the Hamiltonian, we take the same 
operator considered in Ref. [32] for Levy walks on the 
comb. It is 



d 1 

H = -2ihuj [ X— + - 
ox 2 



(13) 



where w is a dimensionless frequency. In the framework 
of the Fokker-Planck equation, it corresponds to the inho- 
mogeneous convection. Here this Hamiltonian describes 
a one-dimensional motion near a hyperbolic point, and it 
has been studied in connection with the Riemann zeros 
and eigenvalue asymptotics Issl-jsol . lioj , scattering of the 
inverted harmonic oscillator [41|, and eigenstates near a 
hyperbolic point (i^ . 

The standard quantum dynamics with the Hamilto- 
nian of Eq. (|13p is linear and coincides with semiclassi- 
cal one. The wave functions and expectation values are 
well behaved values. For example, evolution of the wave 
function is an explicit function of the initial condition 



■^jj{x,t) ^ e 



_ -iHt/h 



(14) 



and defining y — xe '• , one obtains for the second mo- 
ment 



e h 



(15) 



Contrary to this, the semiclassical dynamics of the sys- 
tem described by the FSE (and on the comb, respec- 



dweT^-^Vo (e~"a;) , (16) 



{x\t)) = j y^ro (ye*) ^0 (ye^) dy . (17) 
For example, for the Gaussian initial condition ^0(2;) — 
e"^ /a/tt the expectation value in Eq. (fT7|) diverges at 
t — hir/A. This result was observed in Ref. [44] where 
the quantum motion near the separatrix was studied. 

We can conclude that the fractional time derivative, 
at least for a = 1/2, reflects an effective interaction 
of a quantum system with an additional degree of free- 
dom. As already mentioned, the y direction in the 

comb model was introduced to model the time fractional 

1 

derivative by analogy with the continuous time ran- 

dom walk, where delay times of escapes from the motion 
in the x space are distributed by the power law ~ t~'^/^ 
[Sll . [3^ . [33[ . It is worth noting that in the subdiffusive 
case with a = 1/2 the fractional Fokker-Planck equation 
is identical to the diffusion comb model 32]. In the quan- 
tum case the situation differs essentially from fractional 
diffusion. First of all, the quantum comb model and FSE 
are not identical. As shown here the systems are isospec- 
tral: they have the same singularities for the Green func- 
tions. In the quantum case a "random entrapping" is due 
to the reversible leakage probability of the wave function. 
The decay term in the Green function of Eq. is the 
specific property of the FSE which violates the Hcrmi- 
tian property of Hamiltonian H . For the quantum comb, 
the y space is the environment. In this connection, an 
interesting question arises: how does the time fractional 
derivative, which describes a specific interaction with the 
environment, relate to a possible des crip tion in the frame- 
work of the Lindblad equation [i^, |46| ? Both the quan- 
tum motion on the comb and the FSE introduce new 
nonlinear phenomena in the semiclassical limit, and this 
semiclassical approach differs from one described in the 
framework of the standard Schrodinger equation. 
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